Spin-polaron model: transport properties of EuBg 
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To understand anomalous transport properties of EuBs, we have studied the spin-polaron Hamil- 
tonian incorporating the electron-magnon and electron-phonon interactions. Assuming a strong 
exchange interaction between carriers and the localized spins, the electrical conductivity is calcu- 
lated. The temperature and magnetic field dependences of the resistivity of EuBg are well explained. 
At low temperature, magnons dominate the conduction process, whereas the lattice contribution 
becomes significant at very high temperature due to the scattering with the phonons. Large negative 
magnetoresistance near the ferromagnetic transition is also reproduced as observed in EuBg. 

PACS numbers: 72.10.Di, 71.38.-k, 75.47.Pq 



Hexaboride compounds have been studied extensively 
for their unusual transport and magnetic properties. 
Among those, EuBg has attracted special interest due 
to its exotic magnetic properties, such as two consecu- 
tive phase transitions at low temperature (15.5 K and 
12.6 K), very low carrier density (~ 10^" cm~'^), metal- 
lic ferromagnetism below Tc, and large negative colossal 
magnetoresistance at the transition of 15 K It was 

proposed that the higher transition temperature is a met- 
allization temperature due to an increase in the number 
of itinerant electrons and the lower one is a bulk Curie 
temperature of long-range ferromagnetic (FM) order . 
Sample dependence is an important issue in EuBg as it 
has a very small number of intrinsic carriers. The higher 
transition temperature is seen to be much more sample 
dependent than the lower one 0]- 

Raman scattering measurements 0, Q indicated the 
spontaneous formation of magnetic polarons, involving 
FM clusters of Eu^+ spins just above Tc- It was conjec- 
tured Q that the transition at higher temperature arises 
from the mutual interaction between the magnetic mo- 
ments and the conduction electrons which leads to the 
formation of bound magnetic polarons. The bound mag- 
netic polaron corresponds to a composite object of local- 
ized charge carrier and its induced alignment in a back- 
ground of local moments. On the other hand, Hirsch 
[5i] proposed a model that the magnetism of EuBg is 
driven by the effective mass reduction or the band broad- 
ening upon spin polarization. The origin of this effect 
is the bond-charge Coulomb repulsion which, in a tight 
binding model, corresponds to the 'off-diagonal' nearest- 
neighbor exchange and pair hopping matrix elements of 
the Coulomb interaction. But the parameters for the 
exchange and pair hopping considered in his model are 
unphysical and the Eu 4/-states were treated as delocal- 
ized electrons in contrast to the band structure results 
Pi . The RKKY interaction among the localized 4/ elec- 
trons through itinerant electrons is also considered to be 
an origin of ferromagnetism in EuBgQ- More recently, 
a fluctuation-induced hopping model is proposed as a 
transport mechanism for the spin polaron in a param- 
agnetic background of fluctuating local moments In 
this case, temperature dependence of the resistivity is ob- 



tained to be proportional to 7^5/2 fc^T > Jfm and 
to T for ksT ^ Jfm, where Jfm is coupling between 
the local moments. They claimed that this transport 
mechanism is valid for high-temperature phase of EuBg . 
However, no T^/"^ behavior is observed and moreover the 
crossover temperature is too high to be applicable to 
EuBg. 

The diverse properties of EuBg are expected to come 
primarily from the exchange interactions between the 
carriers and Eu^"*" 4/ local moments. Indeed the large 
energy and field dependence of the spin-flip Raman scat- 
tering peak 3| reflects the substantial FM exchange in- 
teraction {Jcf ^ 0.1 eV) between the carrier and the 
Eu^+ spins in EuBg. The magnetic polaron formation is 
favored by the large ferromagnetic Jcf as a trapped 
charge lowers its energy by polarizing the local moments 
of Eu2+. 

The low-carrier density magnetic systems such as EuBg 
share many properties with the CMR manganites, such 
as the insulator-metal transition concomitant with the 
FM transition, large negative magnetoresistance, and the 
formation of magnetic clusters near Tc- Distinctly from 
manganites, however, Eu-based systems do not manifest 
the structural complexities originating from the strong 
electron-lattice coupling associated with Jahn- Teller ef- 
fect. Due to the structural simplicity, EuBg is an ideal 
system for investigation of the interplay between mag- 
netic and transport properties. An essential aspect for 
a theoretical understanding of the properties of EuBg 
is the temperature dependences of the resistivity. So it 
would be of great interest to study the temperature and 
field dependence of resistivity of EuBg . For this purpose, 
we consider the spin-polaron Hamiltonian taking into ac- 
count the magnetic excitations (magnons in the linear 
spin wave approximation) together with lattice excita- 
tions (phonons): 
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where Cia is the annihilation operator for the conduc- 
tion electron with spin a at site i, and bp are the an- 
nihilation operators for a localized spin (magnon) with 
momentum q and for the phonons with wave vector p, 
respectively, and uj^ and are the magnon and phonon 
frequencies, respectively, ujq ~ JfmS where Jfm be- 
ing the direct FM exchange interaction between Eu 4/ 
spins {S). The last two terms are the interaction between 
the itinerant electrons and localized spins and that be- 
tween the local density of the carrier (rii) and localized 
charge vibrations (phonons), respectively. The non-spin 
flip interaction term {S^{ni^ which just shifts the 

on-site energy is not considered, as we are interested only 
in the transport properties. Without loss of generality, 
the Coulomb interaction can be neglected in the very low 
carrier density limit. 

In EuBg, it was observed that clusters of Eu^+ spins 
are formed via the strong ferromagnetic c-f exchange in- 
teraction 0. Hence one can assume a strong exchange 
coupling which induces the locally ferromagnetically or- 
dered spin clusters in the ground state. Formation of 
magnetic polaron is possible in the low carrier density 
magnetic system if the degrees of spin disorder is suffi- 
cient to localize the carrier, but not too high to prevent 
the local FM alignment As for the electron-phonon 
(e-ph) interaction, the study of Boron isotope effect sug- 
gests that the polarons formed in EuBg are predomi- 
nantly magnetic, with a negligible lattice contribution . 
It is, however, not so certain that the magnetic polaron 
in EuBg do not couple to the lattice degrees of freedom of 
Eu ions. Further, thermal conductivity measurement on 
EuBg reveals that the e-ph scattering is strongly favored 
below Tc to describe the variation of the thermal con- 
ductivity Therefore, to explore the role of lattice in 
the magnetic polaron formation, it is worthwhile to con- 
sider the e-ph interaction in the Hamiltonian. 

To decouple the electron-magnon interaction term, we 
employ the following canonical transformation, 
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^y^-^' (cj^ca + clc,^){a,~ - al.). (2) 



In the presence of e-ph coupling, spread and depth of 
lattice deformation can be studied by using the Lang- 
Firsov (LF) transformation 
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Note that, for weak to intermediate e-ph coupling, the 
variational LF transformation jlll | would give more sat- 
isfactory results than LF transformation. In the present 
work, we consider the LF transformation for simplicity. 
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FIG. 1; The normalized resistivity (arbitrary unit) for EuBe 
in the low-temperature region with a)=15 K, f2=168 K, a=8.0, 
/3=0.5, A=0.004 eV, 6=0.002 eV. symbols represent the 
experimental data in /ificm (Ref. ^3)- 



As a result, the transformed Hamiltonian becomes 
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The magnon part remains unchanged after the trans- 
formations. In the transformed Hamiltonian, only the 
hopping term contains non-linear functions of spin-wave 
operators. As we are interested in the transport prop- 
erties in the very low carrier density limit, we will not 
consider the renormalized interactions between carri- 
ers. To calculate the electrical conductivity, dynamics 
of cosh(2;i — Xj)xjXj and sinh(a;i — Xj)xj Xj in the hop- 
ping term of Eq. Q should be properly treated. At low 
temperature, transport is described by an effective band- 
width with a background of magnons and phonons. The 
effective mass of the carrier increases as a result of the 
interaction with the magnons. The effect of the back- 
ground is included by calculating the thermal average of 
cosh (xi — Xj)xlXj and sinh [xi — Xj)xjXj. 

In the conduction, there are two independent pro- 
cesses: elastic and inelastic. For the elastic conduction 
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FIG. 2: Normalized resistivity vs. temperature T (K) for 
different magnetic fields H (in Tesla). Inset figure shows the 
low-temperature variation of resistivity. 



process, thermal average of the terms mentioned above 
is calculated as in the conventional polaron problem (l2| , 
assuming all magnons and phonons to be independent: 
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and (sinh (x^ — Xj)XlXj) 



-0. HereiV-^s 
the magnon and phonon numbers, respectively, 
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Rj and 



.) denotes the thermal average. 
The temperature dependence in the above thermal aver- 
age comes only through N^^^^ and Np^. If we assume 
at finite temperature that the metallic conduction of the 
elastic process CTg is inversely proportional to the effective 
mass and the effective mass is inversely proportional 
to the effective bandwidth, then 



f7e(r) oc (tcosh(a;, - Xj)x}Xj). 



(10) 



For the inelastic process, the conduction is provided 
by the incoherent hopping with emitting and absorbing 
the magnons and phonons. The conductivity is given by 
the Kubo formula 12], 



(11) 



where n is density of mobile carriers, and w is the tran- 
sition probability rate given by the Fermi golden rule. 
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FIG. 3: Normalized magnetoresistance (MR) vs. tempera- 
ture (T/Tc with T'c=15 K) for different values of the magnetic 
field (in Tesla). The same parameters as in Fig.l are used. 
Inset: MR vs. H (in Tesla) at T=15 K. '-f ' symbols represent 
the experimental data (Ref. 0]). 



with the activation energy A — ji^q'^ql^ql'^ + 
^ -f2p'|J7p'p). For simplicity, if we consider single magnon 
(uj) and phonon (Q) frequency, the total conductivity is 
given by the sum of two independent processes. 



a(T)/a(0) 



exp [-aN{uj) - PN{n)] 
(6V(3AA:BT))exp(-A/A:BT), (13) 



where a = Eq\Vq\^ P = TpWp?, b = etS^n/haiO), 
and 17(0) is the zero temperature conductivity. 

For numerical calculation, we consider lu Tc- As for 
f2, since no inelastic neutron scattering measurements 
have been reported on EuBg, we tentatively use the 
renormalized frequency based on the Einstein phonon 
frequency for LaBg jlj|. If we assume il=168 K for the 
localized mode of Eu^"*" ion, the model calculation repro- 
duces the experimental features. The activation energy 
is a function of u and fl, and also a function of Jcf and 
e-ph coupling. We observed that nature of the resistivity 
is sensitive to the energy parameters A and b. The av- 
erage distortion around Eu site is more than an order of 
magnitude smaller than that in the perovskites |l5j . and 
so we choose /3 to be much less than a. 

Figure 1 provides the calculated resistivity (normal- 
ized to T=0) as a function of temperature. The calcu- 
lated resistivity describes well the qualitative features of 
the experimental resistivity 0, 0, 0, 0. For com- 
parison with the experimental data, we have multiplied 
our results by an arbitrary factor and found very sat- 
isfactory agreement with the experimental results. In 
the present model, the resistivity peak near 15 K sig- 
nifies the crossover from a high temperature insulating 
state with localized and isolated carriers (magnetic po- 
larons) to a low temperature conducting state resulting 
from the overlap of the magnetic polarons. Polaron over- 
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lap is marked by the rapid drop in resistivity at the tran- 
sition. Consideration of only the electron-magnon inter- 
action in the present model can reproduce the resistivity 
peak around 15 K. However, the inclusion of phonons 
within this model makes the agreement better. At very 
low temperature, the contribution of magnons to the con- 
ductivity is more important than that of phonons. With 
increasing temperature, the band conduction dominates 
with enhanced effective mass of the carrier so that the re- 
sistivity increases very rapidly. In contrast, at very high 
temperature, the role of the lattice effect becornes signif- 
icant due to the scattering with the phonons 9], where 
the incoherent hopping dominates. 

Now let us investigate the effect of the external mag- 
netic field on the transport. One can include the effect of 
the external magnetic field H for FM magnons by replac- 
ing the magnon frequency uj^ by (wg + goff/is-ff), where 
QcS is the effective g-factor for the localized spins of Eu^+ . 
<?off(= g* + — takes into account the enhanced effcc- 
five field arising form the exchange interaction between 
the carriers and Eu-4/ local moments in the presence of 
the magnetic field. Here g* is the intrinsic g factor for 
the Eu-4/ spin, x ^nd g are the susceptibility and the g 
factor for the conduction electrons, respectively, and Xc 
is the concentration of conduction electron. Taking the 
temperature independent Pauli paramagnetic suscepti- 
bility for X which is given by the density of states N{Ef) 
of conduction electrons ,2Qi| and Xc = 0.01 per unit cell 
0, one can estimate the value of ges to be nearly 
equal to 6. Then, with the above field dependent Wg-, the 
experimental results of large magnetoresistance can be 
achieved. 

The temperature dependence of the electrical resistiv- 
ity with varying the external magnetic field is presented 
in Fig. 2. The field and temperature dependences of the 



electrical resistivity reveal that the charge transport is 
strongly correlated with the higher transition tempera- 
ture In the presence of magnetic field, the transi- 
tion is broadened and shifted to higher temperature. In 
the inset shown is the resistivity in the low-temperature 
region with varying the external magnetic field rang- 
ing from 0.05-0.5 T. At zero field, the resistance drops 
sharply just above the transition temperature. This re- 
sistivity peak is suppressed by the magnetic field in good 
agreement with the observation 0. 

In Fig. 3, we plot the magnetoresistance {MR—{p{H) — 
p(0))/p(0)) predicted by the present model as a function 
of temperature. With increasing temperature, a large 
negative contribution to MR appears. We obtain the 
negative MR even at temperature higher than Tc which 
is consistent with the experiment [1^1 • Inset of Fig. 3 
provides the normalized MR at T=Tc (15 K) as a function 
of the external magnetic field. A large negative MR near 
the FM transition, as depicted in the inset, is indeed 
observed in EuBg The MR approaches —1 in the 

high field limit, reflecting that the zero field resistivity 
is completely suppressed, which is in excellent agreement 
with experimental features. 

In conclusion, we have proposed the model which ac- 
counts qualitatively well for several anomalous transport 
properties observed in EuBg. The strong exchange inter- 
action (Jc/) and the FM magnons in the present model 
can reproduce the temperature and magnetic field de- 
pendences of the resistivity in EuBg. Large negative MR 
occurs even at temperature higher than Tc, which is con- 
sistent with experiments. 
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